What is claimed is: 

A method for producing an enhanced-resolution image by use of M first low- 
resolution images associated with a scene, each of the M first images being 
represented by a set of uniformly sampled and quantized pixels, a relationship of 
a magnification factor existing between the enhanced-resolution image and the M 
first images, said method comprising the steps of: 

(a) storing the M first images; 

(b) selecting one from the M first images as a prototype image, and referring 
the non-selected (M-l) first images to being as (M-l) second images; 

(c) based on the magnification factor, interpolating extra pixels into the set of 
pixels of the prototype image, and then inferring the values of the extra 
pixels according to the values of neighbor pixels thereof to produce an 
interpolated prototype image, and calculating a respective translation 
between each of the (M-l) second images and the interpolated prototype 
image; 

(d) dividing the translations of the (M-l) second images by the magnification 
factor, respectively, to obtain the modulus with respect to the translations 
of the (M-l) second images, selecting, based on a criteria, one from the 
second images whose related modulus are the same, and referring the 
selected second images together with the rest of the second images whose 
related modulus are not the same to being as N third images, wherein N is 
equal to or less than (M-l); 

(e) down-sampling the interpolated prototype image N times according to the 
respective translation between each of the N third images and the 
interpolated prototype image to produce N fourth images which each 
corresponds to one of the N third images; 

(f) calculating the difference between each of the N third images and the 
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corresponding fourth image thereof; 

(g) adjusting the values of the pixels of the interpolated prototype image 
according to an average of the differences calculated in step (f); and 

(h) repeating step (e) through step (g), until the values of the pixels of the 
interpolated prototype image converge to a satisfactory result, and 
referring the interpolated prototype image whose values of pixels 
converge to the satisfactory result to being as the enhanced-resolution 
image. 

The method of claim 1 , wherein step (e) is performed according to the following 
formula: 

wherein g>> is the Ath forth image at the nth down-sampling process, y(»> is the 
interpolated prototype image been adjusted n times, h is the blurring operator 
defined by a point-spread-function, j is the transformation operator with 
respect to the translation between the Ath third image and the interpolated 
prototype image, is is the down-sampling operator. 

The method of claim 2, wherein steps (f) and (g) are performed according to the 
following formula: 

fw=fw +±f;r t - 1 ((cg 4 - 8t M ) t s)* P ) , 

wherein is the interpolated prototype image to be adjusted at (a?+1) time, 

K is the number of the third images been calculated, g k is the Ath third image, p 
is the de-blurring operator, f~ x is the inverse transformation operator with 
respect to the translation between the Ath third image and the interpolated 
prototype image, |s is the up-sampling operator. 

The method of claim 3, wherein in step (c), the inferring of the values of the 



extra pixels based on one selected from the group consisting of a first order 
interpolation, a third order interpolation and a fifth order interpolation. 

The method of claim 4, wherein in step (c), the translation {LT t ) between the /th 
second image and the interpolated prototype image is calculated according to the 
following steps: 

(c-1) assigning at least one interesting point Pfajt) on the ith second image, 
each interesting point corresponding to a local window (w) and a point 
(w ? v) on the interpolated prototype image; 

(c-2) calculating an absolute difference ( LAD^x.y^u, v) ) between each 
interesting point Pi(xy) and the corresponding point (w,v) thereof 
according to the following formula: 
LAD, (x, y\ u, v) = £ \l { (x + m, y + n) - I Q (u + m, v 4- n)\ 

(c-3) calculating a minimum ( LR t (x,y) ) of the absolute difference 
LAD^x.y.u.v) calculated in step (c-2) according to the following 
formula: 

LR t (x, y) = arg min LAD, (x, y\ u, v) 

(w,v) ; 

(c-4) calculating the translation (£7^) with respect to each interesting point 
Pi(x,y) according to the following formula: 

LT i (x, y) = LRj (x, y) - (x, y) • magnificat ion factor • 

J T 

The method of claim 4, wherein in step (c), the translation ( 0 between the ith 
second image and the interpolated prototype image is calculated according to the 
following steps: 

(c-1) calculating an absolute difference (GAD(u,v)) w ^h respect to the 
whole ith second image according to the following formula: 
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GAD t (u 9 v)= ^\l£x 9 y)-I 0 (u + x 9 v + y) 

wherein the point (w,v) is a point, on the interpolated prototype image, 
corresponding to a start point of the ith second image; 

(c-2) calculating a minimum (GR(i)) of the absolute difference (GAD(u,v)) 
calculated in step (c-1) according to the following formula: 

GR(i) = arg min GAD, (u, v) , 

(»,v) 

and considering Gi?(z) as • 
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